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Abstract. Increasing fidelity is the ultimate challenge of quantum information technology. In addition to 
decoherence and dissipation, fidelity is affected by internal imperfections such as impurities in the system. 
Here we show that the quality of quantum revival, i.e., periodic recurrence in the time evolution, can be 
restored almost completely by coupling the distorted system to an external field obtained from quantum 
optimal control theory. We demonstrate the procedure with wave-packet calculations in both one- and 
two-dimensional quantum wells, and analyze the required physical characteristics of the control field. Our 
results generally show that the inherent dynamics of a quantum system can be idealized at an extremely 
low cost. 
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1 Introduction 

Quantum revival, i.e., periodic recurrence of a wave func- 
tion is a fundamental property of a time-dependent quan- 
tum system. The quality of the quantum revival can be 
determined by a simple overlap between the initial and fi- 
nal states after the ideal revival time, i.e., the revival time 
of the unperturbed system. In this respect, the quality 
corresponds to quantum fidelity [1 . 

Fidelity is the most important measure of the function- 
ality of a quantum information device, and in this respect 
it has been subjected to extensive theoretical work [2 . In 
experiments fidelity is greatly affected not only by even- 
tual decoherence and dissipation due to the coupling with 
the environment, but also by internal imperfections such 
as random irregularities and impurities. The standard ap- 
proach to overcome this problem is to increase the sample 
quality or to make it, by other means, particularly resis- 
tant to imperfections. 

An alternative strategy to improve fidelity, which is the 
topic of this work, is to couple the system to an external 
control field that assists the system to overcome and/or 
to compensate the effects induced by irregularities. Such 
scenarios have attracted significant interest in the design 
of high-fidelity quantum gates in optical lattices [3] . Here 
we show that optimizing the control field with quantum 
optimal control theory [4^5 (OCT) leads to a practically 
complete restoration of the fidelity and thus to the full 
quantum revival even at significantly low energies of the 
control field. We analyze in detail the properties of the 
optimized field and demonstrate the scheme for the time 



evolution of one- and two-dimensional wave packets, re- 
spectively. Our findings have broader applicability to gen- 
eral control problems, where the inherent dynamics can be 
assisted by optimized fields subject to strict constraints in 
terms of the intensity and the frequency range. 



2 Optimal control theory 

Since its formulation ^|5j in the 1980s OCT has increased 
its popularity in chemistry and condensed matter physics [6] . 
The central idea of OCT is to replace trial-and-error type 
learning-loop experiments with a rigorous extension of the 
classical control problem to quantum mechanics. In all 
OCT applications the objective is to find an external time- 
dependent field e{t) that drives the system to the prede- 
fined target, e.g., to a certain quantum state, through the 
solution of the Schrodinger equation. 



.dr. X 



(1) 



where dipole approximation is applied with fi = —y (Hartree 
atomic units used throughout). The central idea is to max- 
imize the target functional 



Ji[V] = (^{r,T)\dW{r,T) 



(2) 



at time t = T. In this study we consider two types of target 
operators: (i) projection operator O = |^^f) (^f| with (^p 
as our target state^ and (ii) a local operator O = Pf{'^) 
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representing our target density t = T. The maximum 
value obtained for Ji is referred to the yield. 

The optimized pulse is exposed to two important con- 
straints. First, the fiuence, i.e., the time-integrated inten- 
sity Fq is kept fixed by applying a functional. 



0.999 



J2\e] = 



-a 



f 

Jo 



dte^{t) -Fo 



(3) 



0.99 



2 



where a is a time-independent Lagrange multiplier. Sec- 
ondly, a spectral filter is applied to cut off the highest 
frequencies, e.g., unrealistic photon energies. This is done 
by multiplying the Fourier-transformed pulse by a filter 
function /(cj), that represents the desired frequency range. 
Thereafter, an inverse Fourier transform of the product is 
taken to obtain the time-resolved filtered pulse (see Ref. [7 
for details). 




x10 



field strength (a.u.) 



Fig. 1. (color online) Yield as a function of the field strength 
Finally, the time-dependent Schrodinger equation [Eq. (P^nd maximum frequency cjmax of the optimized field in a revival 



must be satisfied in the control procedure. This gives us 
yet another functional. 



process of a one-dimensional square quantum well distorted by 
a Gaussian impurity. The initial state is a superposition of the 
two lowest eigenstates in a clean system. 



-2Im / {x{t)\idt 
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mim), (4) 



where the auxihary function x(^) can be regarded as a 
time-dependent Lagrange multiplier [4 . 

Variation of the total functional J = Ji H- J2 + ^3 with 
respect to iZ^, and a leads to the control equations 



idtxit) = H{t)x{t), x(T) = 0^(T), 



1. 
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which are solved iteratively [8l[7] by applying here the 
forward-backward propagation scheme of Werschnik and 
Gross [9j. Typically a converged field is obtained within 
100 . . . 300 OCT iterations. We have applied octopus code 
in all the calculations. 



3 Results 

3.1 Superposition in a one-dimensional well 

We start by considering a single particle with mass m = 1 
in a one-dimensional square quantum well with infinite 
boundaries and a length L = 20. The energy eigenstates 
are given by En = 7r^n^/(2mL^) = (7r^/800)n^ with n = 
1, 2, . . . By expanding the wave function and equalizing 
the phase factors it is straightforward to show pT| [T2 i rr3 ] 
that the exact quantal revival time for any wave function 
IS given by Trev = 4.mL^/7T = IGOO/tt. However, if the 
wave function is a superposition of two eigenstates with 
En and Em^ the first full revival occurs already at = 
27r/ {Em — En). 



In our first example our initial state is a superposi- 
tion of the lowest two eigenstates, so that the first revival 
time with |(lZ^(x,Trev)|^(x,0))|' = 1 is T^l = Trev/3 = 
1600/(37r). Starting from the same state but supplying 
our quantum well with a Gaussian impurity Vimp{x) = 
/3exp[— (x — xo)^/7^] with /3 = 0.1, 7 = 1, and xq = 
2.5 leads to a reduced overlap, \{^\x^Trev)\^{x^O))f = 
95.9%. Here and I^Z^") correspond to distorted and 
clean systems, respectively. 

Our task now is to restore the full quantum revival by 
optimizing a control field with a target = ^{x^O) = 
lZ^(x,T). Figure [1] shows the obtained yield (overlap) as 
a function of the maximum allowed frequency cjmax and 
the field strength of the initial constant field e (so that 
the fiuence Fq = T^ev^'^ is kept fixed in the OCT proce- 
dure). Generally, Fig. [T] demonstrates that the optimiza- 
tion works in the desired manner, so that we find a signif- 
[mcant increase in the yield even up to > 99.9%. The price 
to pay is embedded in the allowed intensity and frequency 
of the field as analyzed in the following. 

As shown in Fig.[T]the yield is strongly increased when 
^max is increased above ~ 0.15. It can be expected that 
the " critical" frequency is related with the energy gaps in 
the spectrum. Now the critical frequency is clearly above 
the first few energy gaps, so it seems likely that (de- 
) excitations at relatively high levels are needed in the dis- 
torted system to produce a high overlap with the initial 
state. The coupling with the optimized field with the en- 
ergy levels is studied in more detail in Sec. 13.21 We point 
out that at small cutoff frequencies (cjmax ^0.15) the yield 
is smaller than in the distorted system without a control 
field; with a limited frequency range, and especially when 
supplied with a high field strength, the field is causing 
harm in the system despite the optimization. 

In a reasonable frequency range, the required fiuence 
to reach a high overlap is very small. For example, 99% 
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yield can be obtained with Fq ^ 2 x 10~^, which is only 
about 5% of the first excitation energy Ei — Eq = 0.037. 
The energy carried by the field is also much smaller than 
the difference in the energy levels of the system with and 
without the impurity. Furthermore, in conventional con- 
trol problems where, e.g., charge transfer or a desired ex- 
citation is optimized, the typical intensities are larger by 
several orders of magnitude despite similar system char- 
acteristics p!4[[T5) [T6] . Overall, we may conclude that in 
processes where the inherent dynamics gives a high initial 
overlap (as in the partial revival), OCT finds the route to 
almost 100% yield with an extremely low price. 

3.2 Gaussian packet in a one-dimensional well 

Next we study the revival of a Gaussian wave packet in 
the same one-dimensional quantum well. The wave packet 
is defined as ^{x^O) = exp(— (x + xo)'^)/^/7^ with xq = 
— 3V^, and the impurity is the same as before [see the 
inset of Fig. [2fd)]. As the wave packet consists of a large 
ensemble of eigenstates the first revival will occur at the 
universal revival time T^ev = 1600/7r as discussed in the 
previous section. This is set as the time of the control 
field, and the target is now the (single-particle) density 

PF(^,T = Trev) = 1^^(^,0)12. 

Figure [2] summarizes the OCT results obtained with a 
threshold frequency cjmax = 1 and fluence Fq = 0.00203. 
Figure [2](a) shows the overlap with the initial wave packet 
as a function of time in the case of (i) a clean system (solid 
line), (ii) a distorted system with the impurity (dotted 
line), and (iii) a distorted but controlled system (dashed 
line). During the last stages of the propagation [Fig. [2fc)] 
the controlled wave packet closely follows the ideal system 
and reaches 99.9% overlap, whereas the uncontrolled evo- 
lution leads to about 88%. Thus, the optimization seems 
to find the correct "path" in time, so that the field contin- 
uously corrects the evolution toward the ideal one, even 
though the target functional is not time-dependent. 

The optimized control field and its Fourier spectrum 
are shown in Figs.[2fb) and (d), respectively. The strongest 
two peaks in the spectrum at = 0.025 . . . 0.05 correspond 
to the transition between the first two levels with AE = 
0.037, and the peak at ~ 0.15 is likely to correspond to 
the transition between the 5th and 6th states with AE = 
0.13. Figure [2fe) shows the projections of the propagated 
single-electron orbitals (of the system with an impurity 
and a control field) to the initial orbitals (of the clean 
system). Here the excitation processes between the 1st and 
the 2nd, and, on the other hand, between the 5th and the 
6th states are clearly visible. At the end of the pulse the 
initial occupations are reached as the initial wave packet 
is reconstructed in the controlled revival process. 

In Fig. [3] we show the effect of the field strength and the 
maximum allowed frequency on the obtained yield in the 
controlled revival of a Gaussian wave packet. The criti- 
cal field strength (of the initial constant field) to reach 
a yield close to 0.999 - after the optimization within a 
fixed fiuence - is e ~ 1.2 x 10~^ corresponding to Fq ~ 
0.00073. As a function of cjmax we find two steps at around 




0.5 1 1.5 0.01 0.1 1 
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Fig. 3. Yield in the controlled revival of a one-dimensional 
Gaussian wave packet as a function of (a) the field strength 
(without a filter, i.e., cjmax oo) and (b) the maximum al- 
lowed frequency when the field strength of the initial pulse is 
e = 0.002. 

0.004. . .0.005 and 0.15, respectively. They correspond to 
the excitations 1^2 and 5^6 and confirm the impor- 
tance of those processes in the controlled evolution. 

Next we consider the effect of the impurity characteris- 
tics on the controllability. It is expected that a narrow im- 
purity as the one in our example [inset of Fig.[2jd)] probes 
higher excitations than a broader impurity, and hence the 
dynamics and control procedure are qualitatively different 
in those cases. To test this, we increase the width param- 
eter 7 of the Gaussian impurity from 1 to 3 in Vimp{x) 
(see Sec. 13. ip and decrease the height P from 0.3 to 0.2. 
This change reduces the quality of the revival from 88% to 
80%. However, there is no significant change in the control 
process: the yields are generally comparable, and the step 
in Fig.[3fa) is at the same place. The yield as a function of 
the threshold frequency is also similar, although the step 
structure is not so prominent. 



3.3 Gaussian packet in a two-dimensional well 

Finally we move our attention to a two-dimensional ex- 
ample. This extends the physical relevance of the pro- 
posed scheme to two-dimensional electron gases having 
a variety of applications in, e.g., quantum Hall and quan- 
tum dot physics [17 . We consider a square quantum well 
with hard- wall boundaries and side length L = 12 [see 
Fig. m^a)]. The initial Gaussian wave packet is given by 
^{x,y,0) = 5 exp(-(x + xq)^ + (^ + yo)^)/{25^) with 
(5 = 0.7, xo = 1, and yo = 2. We note that due to symme- 
try the wave-packet propagation is separable to x and y 
components. However, we expose the system to ten ran- 
dom repulsive impurities visualized in Fig. ID^b) (two of 
them merged together), so that the problem is no longer 
separable. 

Time-evolution of the wave packet in the distorted po- 
tential leads to 71% revival quality. The overlap as a func- 
tion of time is shown in Fig. [5ja). To correct the process 
we add a control field with two independent components 
polarized in x and y directions, respectively. As in the 
previous example, the target is the density of the initial 
(single-electron) wave packet at t = T = Trev The total 
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Fig. 2. (color online) Controlled revival of a one-dimensional wave packet in a square quantum well, (a) Time-dependent overlap 
of the propagated wave packet with the initial one in a clean system (solid lines), in a distorted system (dotted lines), and in 
a distorted by controlled system (dashed lines). In (c) the last part of the propagation is zoomed, (b) Optimized control field. 

(d) Fourier spectrum of the control field. The inset shows the external potential with the initial wave packet and the impurity. 

(e) Projections of the propagated single-electron orbit als to the initial ones. The maximum allowed frequency is ojmax = 1 and 
the fluence is fixed to Fq = 0.00203. 



field fluence is fixed to Fq = 0.07348 and the maximum 
frequency allowed in the optimization is cjmax = 0.5. With 
these constraints OCT produces a control field shown in 
Fig. [5jb) that drives the system to 97% overlap. Thus, 
the control procedure works well also in a two-dimensional 
system. 

Figures [Hc-e) show the single-electron densities dur- 
ing the time evolution at a half of the revival time T/2, at 
3T/4, and at the revival time T. Interestingly, at T/2 the 
distorted packet is rather close to the ideal one as shown as 
a high overlap in Fig. [5fa). Towards the end of the evolu- 
tion, however, the controlled procedure starts to resemble 
the ideal one, whereas the distorted process becomes very 
different. Hence, it seems that the essential adjustments 
in the controlled evolution always occur during the last 
stages of the process. The relevance of the earlier stages 
in the "system preparation" could be studied by consider- 
ing partial control fields affecting the system only during. 



e.g., the last 10% of the time evolution. In addition, it 
would be worthwhile to examine if the ideal evolution (in 
a clean system) can be followed continuously in time by 
applying a time-dependent target functional [7 . These is- 
sues are beyond the scope of the present study and subject 
to future research. 

Finally we comment on the simplification of using a 
fixed and static impurity configurations in the above ex- 
amples. The obtained control field is found only for a spe- 
cific configuration; a general field able to control a large 
ensemble of potentials would be impossible to find in prac- 
tice. Nevertheless, we point out there are several physical 
situations, e.g., in semiconductor heterostructures, where 
the impurity configuration in a given experimental setup 
can be found. For example, in a quantum-dot experiment 
in Ref. [18j the size and position of a migrated impurity 
ion was accurately determined in a system that was oth- 
erwise close to an ideal two-dimensional harmonic oscil- 
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Fig. 4. (color online) (a-b) External potentials of clean and dis- 
torted square quantum wells, respectively, (c-e) Single-electron 
densities at times t = T/2, 3T/4, and T of the propagation in 
a clean, distorted, and controlled system, respectively. 
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Fig. 5. (color online) (a) Overlap as a function of time in a 
revival process of a two-dimensional Gaussian wave packet in 
a square quantum well, (b) Optimized two-component (x and 
y) control field leading to 97% overlap. Here Fq = 0.07348 and 

^^max = 0.5. 



lator. In such a situation, the Hamiltonian is known and 
optimization of the field to obtain a predefined target is 
possible. 



4 Summary 

To summarize, we have analyzed quantum revival pro- 
cesses of single-particle states and Gaussian wave pack- 
ets in one- and two-dimensional quantum wells. We have 
shown that the quality of the quantum revival in a realis- 
tic (distorted) system can be greatly improved by coupling 
the system to an external control field. The control field 
can be optimized with quantum optimal control theory 
to maximize the overlap between the initial wave func- 
tion and the time-propagated one at the revival time - 
within predefined constraints on the fluence and maxi- 
mum allowed frequency. We have analyzed how the field 
constraints affect the obtained yields, and found that very 
low intensities are sufficient to obtain yields above 99% 
in one-dimensional systems. The threshold frequencies to 
achieve high yields can be associated with excitation en- 
ergies in the system. Our procedure has broader implica- 
tions to general control problems, where the objective is to 
supplement the inherent dynamics disrupted by irregular- 
ities. The demonstrated applicability to two dimensions 
shows that the proposed approach could be used to in- 
crease quantum fidelity in, e.g., quantum Hall devices. 
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